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ABSTRACT 

Since the photons that stream from QSOs alter the ionization state of the gas they traverse, any 
changes to a QSO's luminosity will produce outward-propagating ionization gradients in the surrounding 
intergalactic gas. This paper shows that at redshift z ~ 3 the gradients will alter the gas's Lyman-a 
absorption opacity enough to produce a detectable signature in the spectra of faint background galaxies. 
By obtaining noisy (S:N~ 4) low-resolution (~ 7 A) spectra of a several dozen background galaxies in a 
R ~ 20' field surrounding an isotropically radiating 18th magnitude QSO at z = 3, it should be possible 
to detect any order-of-magnitude changes to the QSO's luminosity over the previous 50-100 Myr and 
to measure the time tQ since the onset of the QSO's current luminous outburst with an accuracy of 
~ 5 Myr for tQ £ 50 Myr. Smaller fields-of-view are acceptable for shorter QSO lifetimes. The major 
uncertainty, aside from cosmic variance, will be the shape and orientation of the QSO's ionization cone. 
This can be determined from the data if the number of background sources is increased by a factor of a 
few. The method will then provide a direct test of unification models for AGN. 

Subject headings: quasars: general — galaxies: high-redshift — intergalactic medium — black hole 
physics 



1. INTRODUCTION 

The 10 8 -M Q black holes that lie at the heart of nearby 
bulge galaxies are believed to have accreted much of their 
mass in their youths when they shone briefly and brightly 
as QSOs (e.g., Yu & Tremaine 2002). A great deal could be 
learned about the physical processes that produce super- 
massive black holes if we could observe how their bright- 
nesses varied during this time. For example, a black 
hole that radiates at a fraction I = L/Le of its Edding- 
ton luminosity and accretes mass with radiative efficiency 
e = L/Mc 2 requires t s ~ M BH /M BH ~4x 10 7 ;(e/0.1) yrs 
to change its mass significantly. Showing that QSO out- 
bursts have a much shorter duration would confirm the 
popular belief (e.g., Haehnelt & Kauffmann 2000; Cava- 
liere & Vittorini 2000; Wyithe & Loeb 2002; Di Matteo 
et al. 2003) that super-massive black holes build up their 
masses through numerous accretion episodes. Various the- 
oretical attempts to explain the observed correlation be- 
tween a black hole's mass and its bulge's velocity disper- 
sion (e.g., Silk & Rees 1998; Adams, Graff, & Richstone 
2001; Burkert & Silk 2001; Miralda-Escude & Kollmeier 
2004) each postulate different physical mechanisms that 
quench the QSO's luminosity when it reaches a certain 
level. To the extent that these mechanisms operate on 
different time scales, measuring the duration of QSO out- 
bursts should help us distinguish between them. 

The challenge is to figure out from a few short nights or 
decades of observations how a QSO's brightness changed 
over the preceding few million years. Investigators have 
resorted to a number of indirect schemes. By making as- 
sumptions about the link between QSOs and either local 
black holes or massive potential wells at the redshift of 
observation, several authors have deduced that individ- 
ual QSOs are unlikely to have shone for fewer than lMyr 
or more than lOOMyr (e.g., Haehnelt, Natarajan, & Rees 
1998; Richstone et al. 1998; Martini & Weinberg 2001; 
Hosokawa 2002). Unfortunately the uncertain underlying 



assumptions do not seem to allow a much more precise 
limit on the typical QSO lifetime. In any case statistical 
arguments like these will never be able to tell whether a 
QSO's total radiative lifetime of (say) 10 7 yr occurred in 
one contiguous chunk or was instead split among numer- 
ous shorter bursts, though this distinction could be crucial 
in developing a physical picture of black-hole formation. 

This paper describes a method that is somewhat less 
indirect and that can provide a rough indication of how 
any given QSO's luminosity varied over the 50-100 Myr 
preceding the time of observation. The simple idea behind 
the method is illustrated in Figure 1. Although we cannot 
detect the photons a QSO emitted in the past, we can de- 
tect their effect on its surroundings. As hydrogen-ionizing 
photons from a QSO propagate outwards, they destroy 
neutral hydrogen in the intergalactic medium and reduce 
the number of Lyman-a absorption lines in the spectra of 
background objects. The change in the Lyman-a opacity 
at radius r should therefore provide some indication of the 
QSO's luminosity at the earlier time t ~ r/c. The rest of 
the paper works out this simple idea in more detail. § 2 
presents a brief review of the relevant intergalactic physics 
and justifies two assumptions that will be needed later. 
§§ 3 and 4 work out the effect of changes in the QSO's 
luminosity on background galaxies' spectra. The next two 
sections discuss the significance with which the effect can 
be detected: § 5 treats the uncertainties from a theoretical 
point-of-view and concludes that cosmic variance will be 
the primary problem, while § 6 discusses the spectroscopic 
exposure times that will be necessary on a 10m telescope 
and argues that neither continuum removal nor interstel- 
lar absorption lines will be a major obstacle. The results 
from the preceding sections are brought together in one 
place in § 7, which presents a sample analysis of a simu- 
lated QSO with a known radiative history. Some readers 
may wish to skip directly to this section. § 8 discusses the 
time resolution that can be achieved with this technique. 
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My main conclusions are reviewed and criticized in § 9. 

I should state at the outset that this is not the only 
method for constraining the radiative histories of individ- 
ual QSOs. Readers may judge for themselves the relative 
merits of the alternatives that are listed (e.g.) in the ex- 
cellent review by Martini (2003). Nor is the idea behind 
the method new. Jakobsen ct al. (2003) have used it, for 
example, to estimate an age of t > 10 7 yr for the QSO 
Q03022-0023 from the lack of absorption lines in the spec- 
trum of a neighboring QSO. Schirber, Miralda-Escude, & 
McDonald (2004) and Croft (2004) applied a similar anal- 
ysis to a larger sample of QSO pairs and inferred signif- 
icantly shorter lifetimes. What is new, as far as I know, 
is the assumption that the background sources will be nu- 
merous faint galaxies rather than a single bright QSO. This 
complicates the analysis in a number of ways but yields a 
considerably more detailed view of the foreground QSO's 
radiative history. 




Fig. 1. — Schematic view of the underlying concept. QSOs emit 
H-ionizing radiation that will destroy some of the HI in the inter- 
galactic medium out to a maximum radius ctQ , where tQ is the time 
since the onset of QSO activity. HI in the intergalactic medium pro- 
duces absorption lines in the spectra of background objects. The 
absorption is represented by a broken line between the objects and 
the observer; breaks in the line mark positions with significant HI 
absorption. By obtaining spectra of numerous galaxies that lie be- 
hind a high-rcdshift QSO, one can (a) map the spatial distribution 
of intergalactic HI around the QSO, (b) use the shape of the region 
with lowered HI absorption to estimate the rough geometry of the in- 
tergalactic volume that has been hit by the QSO's radiation, and (c) 
deduce the QSO's lifetime and the level of anisotropy in its emission. 



2. PRELIMINARIES 

The low observed level of Lyman-a absorption from in- 
tergalactic gas at z ~ 3 implies that the gas must be al- 
most completely ionized, with perhaps only one neutral H 
atom per million (Gunn & Peterson 1965). As a result 
the typical hydrogen ionizing photon will travel far before 
it is absorbed, ~ 50 proper Mpc (Am = 0.3, Oa = 0.7, 
h = 0.65; Madau, Haardt, & Rees 1999), and one can 
safely assume that intergalactic gas is optically thin on 
length scales significantly smaller. Radiation from distant 
sources will permeate this optically thin gas at a roughly 
uniform level Jb g , ionizing residual hydrogen atoms at a 
rate T^nni where hhi is the neutral hydrogen density 
and 

/>oo 

T* = / dE4nJ hg {E)ai{E)/E (1) 

J13.6cV 

is an integral over energy of the photon number density 
times the hydrogen photoionization cross-section cjj. For 
plausible intensities of the background radiation field, e.g., 

J b gH ~ 10" 21 ' 3 (/i^/13.6eV)" 1 - 8 ergs- 1 cm- 2 sterad" 1 Hz" 1 , 

(2) 



and for densities near the cosmic mean, photon absorption 
will be the dominant ionization pathway for hydrogen, and 
the neutral fraction r\ of intergalactic gas will adjust itself 
until the photoionization rate is equal to the recombination 
rate aHiin H nn e : 

V= T — — p ■ (3) 

nm + nmi 1 7 

The recombination coefficient has a temperature depen- 
dence that is well fit by the expression ami — 2.11 x 
10-14^-0.7^ +T o.7)-i cm 3 s -i where Tg ig the tempera _ 

turc in units of 10 6 K (Cen 1992). 

Now consider what happens to intergalactic gas when 
it is hit by a blast of ionizing radiation from a nearby 
QSO. The photoionization rate increases by an amount 
T^, given by equation 2 with the QSO's radiation field 
Jq replacing Jbg, and the neutral fraction rj falls to its 
new equilibrium value ?7 now = ?7 old r^ s /(r^ g + T®) on the 
time scale \nmjnm\ ~ or ~ 10 4 yr if Jq is equal 

to the background intensity Jbg- After the blast subsides, 
recombination will raise the neutral density back to its 
previous level on the same time-scale. 

In contrast to the potentially large swings in the gas's 
neutral fraction, any changes to its temperature should 
be imperceptibly slight. Although the gas will warm as 
photo-electrons collide with other particles and distribute 
their kinetic energy, the change in its total thermal energy 
will be negligible: photo-electrons are necessarily as rare 
as neutral atoms (~ lppm) and their typical kinetic energy 
at ejection 

em = (r?)- 1 / dEATTj Q (E)a l (E)(E-E th )/E(A) 

JE th 

= E th /(2 + s) ~ 3.6cV for J Q (E) cx E~ s , 
s = 1.8, ai(E) cx E~ 3 , E th = 13.6eV 

is not very different from the energy per particle in the 
~ 20000K undisturbed gas. The gas received its energy of 
about an eV per particle when it was almost completely 
ionized at earlier times, and its temperature will hardly be 
affected by giving a few eV to the particle per million that 
remained neutral. 2 This shows that the temperature of the 
gas will not rise appreciably while its ionization balance 
adjusts to the increased ionizing background. Afterwards 
the gas will be heated by photoionizations at the same 
rate as before — the increase in photoionization rate per 
HI atom will be exactly compensated by a decrease in the 
density of HI atoms — and so the equilibrium temperature 
will be the same in parts of the IGM that are and are not 
illuminated by the QSO's radiation. 

2 I am neglecting the possibility that the gas around the QSO 
may be heated by photo-ionization of atoms other than hydrogen. 
In order to cause a significant temperature change, the product of 
the photo-ionized atom's number density tij and typical ejection 
energy ex must be comparable to nukT. Hell will satisfy this con- 
dition before it is reionized; it is abundant and its photo-electrons' 
ejection energies are large due both to its high ionization potential 
(sec equation 5) and to its optical thickness. The latter means that 
essentially all photons more energetic than the ionization threshold 
will be absorbed, not merely the lower energy photons whose ioniza- 
tion cross-section is greatest (Abel & Hachnclt 1999). By rcdshift 
z ~ 3 the reionization of Hell should be nearly complete, and my ne- 
glect of this heating source is justified. It would not be if the redshift 
were much greater. 
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Taken together, the results reviewed in this section jus- 
tify two assumptions that I will adopt for the remainder 
of the paper, (a) If a QSO's luminosity towards solid an- 
gle at time t is L(Q, t), then the intensity of the QSO's 
radiation field at position r, £1 at time t + |r|/c will be 
proportional to L(fl,t)/r 2 . In other words, intcrgalactic 
gas at larger distances will not be significantly shielded 
from the QSO's radiation by intergalactic gas at smaller 
distances, (b) The intergalactic temperature will not be 
affected by changes in the intensity of ionizing radiation. 
This implies first that a blast of radiation from a QSO 
will alter the ionization balance of intergalactic gas but 
not its spatial or velocity structure, and second that any 
changes in neutral fraction will be precisely proportional 
to the change in the ionizing radiation density (equation 3 
with constant otHii) averaged over the last ~ 10 yr. 

3. MEAN TRANSMISSIVITY VERSUS J v 

We will be able to measure changes in a QSO's ioniz- 
ing history with this method only if we can detect spatial 
changes in the surrounding density of neutral hydrogen. 
Intergalactic HI density is normally measured by fitting 
Voigt profiles to the numerous Lyman-a absorption lines in 
the spectra of background QSOs. This approach is not fea- 
sible when the background sources are faint galaxies, since 
individual Lyman-a forest absorption lines are hopelessly 
blended and confused in their noisy low-resolution spec- 



tra. Instead we can only hope to measure / = 



the mean transmissivity along a spectral segment whose 
length of a few A is many times larger than the typical 
absorption-line spacing but comparable to the instrumen- 
tal resolution. Although a reliance on / is forced on us, 
/ has two advantages over rim as a probe of the ionizing 
background radiation. First, it receives significant con- 
tributions from the parts of the IGM with middling HI 
optical depths TL yct ~ 1 whose response to changes in the 
ionizing radiation are easiest to measure. In contrast nni 
is dominated by systems with large optical depths, and 
since e~ T ~ e~ 2r ~ for large r, big changes in the col- 
umn densities of optically thick systems can be hard to 
detect. Second, the value of nni along any particular line- 
of-sight is strongly affected by whether the line happens to 
pierce an especially dense system. This adds to fini large 
random fluctuations that may obscure underlying changes 
in the radiation field. / weights more evenly across sys- 
tems of different column densities and is consequently less 
affected by chance fluctuations in the density of matter, a 
point that we will elaborate upon in § 5 below. 

Let us work out, then, how / responds to changes in 
the ionizing background. As argued in § 2, r at fixed to- 
tal hydrogen density is inversely proportional to the ioniz- 
ing radiation intensity J, and so if the radiation intensity 
changes from J to b J, with b an arbitrary constant, the 
new mean transmissivity will be 



= / 

Jo 



drP(r)e 



-r/b 



(5) 



where P(r) is the Lyman-a optical-depth distribution that 
existed before the change. 

To estimate P{t), I converted published Lyman-a Voigt 
profile lists for 5 QSOs 3 near z = 3 into lists of the ob- 

3 HS1946+7658 from Kirkman & Tytler (1997; z = 2.994) and 



served Lyman-a optical depth r vs. comoving distance, 
applied a QSO-dependcnt scaling to every r to make each 
QSO's Lyman-a forest have the same mean transmissiv- 
ity / — 0-67 appropriate to z = 3.00 (McDonald et al. 
2001), and finally constructed a histogram of the resulting 
t values. 

The top panel of figure 2 shows the dependence of the 
mean transmissivity on the radiation intensity, calculated 
numerically from equation 5 using a spline fit to this 
histogram as an approximation to P(t). Also shown is 
Af = (df/d In Jq), with Jq = J v — Jb g the QSO's con- 
tribution to the radiation _field. This provides some in- 
dication of how strongly / responds to large changes in 
Jq. The bottom middle panel presents the </„-/ relation- 
ship in a slightly different way, as the mean transmissivity 
as a function of distance r from an isotropically radiating 
z = 3.0 QSO with AB magnitude m 9 i 2 = 18 or m 9 i 2 = 20 
at rest-frame 912A that has been radiating forever at a 
constant rate. This panel assumes that the background 
radiation field Jb g (equation 2) would have been spatially 
uniform in the absence of the QSO, which implies that the 
actual radiation field is 



J(r) 



1 



where 



' eq 



d L (z) 



fv/Jhg 



Jbg 



4tt(1 + z) 



1/2 



(6) 



(7) 



/„ = io-°- 4 ( m 9i2+48.60) is thc flux from the qso received 
on earth at wavelength (1 + z)x912A and d^z) is the 
QSO's luminosity distance. 

The top panel implies (for example) that if the mean 
transmissivity in a given region can be estimated with 
an uncertainty of Uf ~ 0.04, then a significant change 
to the QSO's luminosity (Aln Jq ~ 1) will be marginally 
detectable if the region's distance to the QSO implies an 
expected mean transmissivity 0.72 ?S f ?S 0.93. The bot- 
tom panel shows that the required distance is 12 ;$ r 2 
proper Mpc for a QSO with ragi2 = 18. Adopting the 
crude approximation tdoiay ~ — 2r/c leads to the prelim- 
inary guess that the method should provide reasonable 
constraints on the QSO's luminosity for time delays of 
— 10 £ t £ —80 Myr. A more careful treatment is deferred 
until § 7. 

4. TIME-DELAY SURFACE 

The actual situation is slightly more complicated than 
Figure 1 suggests, since light from the background sources 
does not pass the QSO instantaneously. Instead it encoun- 
ters material that lies behind the QSO before it encounters 
material that lies in front, and as a result the observed in- 
tergalactic absorption from material behind the QSO will 
be sensitive to the QSO's luminosity at an earlier time. If 
we define t = as the time when the QSO emitted the 
light that is just now reaching earth, and if we place thc 
QSO at thc origin of a polar coordinate system where R 
measures proper displacements along the plane of the sky 
and z measures proper displacements in the redshift di- 
rection, then the light from a background galaxy passes 
through the intergalactic gas that lies a proper distance z 

Q0636+680, Q0956+122, Q0302-003, and Q0014+813 from Hu et 
al. (1995; z = 3.180, 3.288, 3.294, 3.366) 
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Fig. 2. — Top panel: the sojid line shows dependence of intcrgalac- 
tic Lyman-a transmissivity / = (e~ TL y Q ) at z = 3 on the intensity 
J v of the ionizing radiation field. The transmissivity / ~ 0.67 is 
appropriate to the actual background ionizing-radiation field Jj-, g . 
Increases in the radiation intensity destroy neutral hydrogen and 
increase the transmissivity. The curve shown was calculated with 
equation 5 from Voigt-profilc fits to the Lyman-a forest at z ~ 3 
(see text). The dashed line shows A/ = (df /djy)^^ — Jbg)> which 
is roughly the change in / that would result if the perturbation to 
the radiation field (i.e., the difference between J„ and Jb g ) were dou- 
bled. Bottom panel: the expected intcrgalactic transmissivity as a 
function of distance from a constantly shining QSO with apparent 
AB magnitude mgi2 = 18 (solid line) or mgi2 = 20 (dotted line) at 
z = 3. The QSO's radiation adds to the background radiation Jb g , 
which we are taking to be J^g = 5 X 10~ 22 erg s~ 
at 912A, and increases the transmissivity in its vicinity. 



behind the QSO at the time t = —z/c, and the ionization 
balance of the gas at this time is sensitive to the QSO's 
luminosity at the earlier time 



t I (R,z)=[-z-(R 2 + z 2 ) 1 / 2 ]/c. 



(8) 



Figure 3 shows the parabolic contours of this function in 
astronomically useful units. The figure implies, for ex- 
ample, that the ionization balance of the intergalactic gas 
that lies 5 Mpc behind the QSO and 2 Mpc to the left will 
reflect the QSO's luminosity at time t = —34 Myr, while 
the material that lies directly in front of the QSO will have 
an ionization balance that reflects the QSO's luminosity at 
t — (i.e., its observed luminosity). 

Now the field of view of a large optical imagers is around 
40', or - 19.9^ proper Mpc at z = 3 (fi M = 0.3, 
Q\ = 0.7), and the largest optical multi-object spectro- 
graph (Dressier, Sutin, & Bigelow 2003) on an 8m-class 
telescope is not much smaller. This makes it relatively 
easy to obtain spectra of high-rcdshift galaxies throughout 
anr~ 20' region centered on a QSO at z = 3, probing 
the ionization balance of the IGM in the region of this 
plot with -10 < R < 10 Mpc. Because more than 1000 
background galaxies with redshift z £ 3.0 and magnitude 



Fig. 3. — The time-delay surface. A QSO at the origin 
emits hydrogen-ionizing radiation into the surrounding intergalactic 
medium. Concentric circles show the distance that has been trav- 
eled traveled by the photons emitted by the QSO at the earlier times 
t = — 10, -20, -30, -40, -50 Myr. Photons from background sources at 
z — ► oo traverse this region from top to bottom as they travel to earth 
at z — ► — oo. At each point along the way they pass photons that 
were emitted by the QSO at a different time in its past. The shaded 
contours show the emission time of the QSO photons that were il- 
luminating a given intergalactic region when the observed photons 
from a background source passed by. The strength of the Lyman-a 
absorption lines from this region in the spectrum of the background 
source will depend on the QSO's luminosity at the indicated time. 
Averaging together the absorption lines from large regions with the 
same time delay (e.g., the region marked A) can provide a constraint 
on the QSO's luminosity at that time. Comparing many such regions 
allows one to chart historical fluctuations in the QSO's luminosity. 
A particularly promising strategy is to estimate the luminosity at 
the time corresponding to region A (~ —90 Myr) by comparing the 
absorption lines in A to the absorption lines in a similar region (B) 
on the opposite side of the QSO that is illuminated by QSO photons 
emitted at t ~ 0. Most systematics will cancel out in this binary 
comparison. See text. 



TZ < 25.5 will be found in a field of this size with the 
"UV drop-out" technique (e.g., Steidel et al. 2003), the 
Lyman-a absorption lines in these galaxies' spectra can 
in principle provide an enormously detailed view of the 
IGM's ionization balance near the QSO. 

Assume, then, that we have measured the intergalactic 
Lyman-a transmissivity / = exp(— TL yct ) throughout the 
region of the figure with — 10 < R < 10 Mpc. A number of 
ways to estimate the evolution of the QSO's ionizing lumi- 
nosity suggest themselves immediately. One is to divide 
the region behind the QSO (i.e., the region with z > 
in figure 3) into a number of bins whose edges align with 
contours of the time delay surface, then see how the mean 
transmissivity / bs in each compares to / exp , the expected 
transmissivity if L(t) were constant. Region A in figure 3 
is one such bin. The problem is estimating / eX p- This 
suggests a slight variation: compare / h s in each bin to 
/ con , the mean transmissivity in a bin of the same shape 
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located on the opposite side of the QSO. See (e.g.) regions 
A and B in the figure. Since A is illuminated by photons 
emitted at t ~ —80 Myr and B is illuminated by pho- 
tons emitted at — lOMyr £ t < 0, the difference in their 
mean transmissivities will be sensitive to the any differ- 
ences in the QSO's luminosity at those two times. This 
approach is also not ideal, since the mean transmissivity 
in the control region (B) is an unnecessarily noisy indica- 
tor of the QSO's luminosity at time — lOMyr < t < 0; only 
a small fraction of the volume illuminated by the QSO's 
luminosity at — lOMyr £ t < falls in region B. The best 
approach is probably to find a maximum-likelihood fit of 
the data to an appropriate surface. Since the results of 
unbinncd maximum-likelihood fits are harder to present 
in a simple graphical way however, I will continue with a 
paired-bin analysis for the remainder of this paper but will 
make one assumption that should cause the estimated un- 
certainties to more closely resemble those of a maximum- 
likelihood fit: I will assume that the uncertainties in mean 
transmissivities of the control bins are negligibly small. 
In other words, I will assume that the uncertainty in the 
mean transmissivity in region B of figure 3 is small com- 
pared to the uncertainty for region A. The justification is 
that random fluctuations in B can be removed to a large 
extent by analyzing the large volume illuminated by the 
QSO's luminosity at -lOMyr < t < 0. 4 

The fact that so large a region is illuminated by light 
emitted at t ~ is crucial to the success of this approach, 
since it allows historical changes in the ionizing luminosity 
to be measured by comparing regions that are distributed 
symmetrically about the QSO. Various systematics (e.g., 
the r -2 decrease in flux from the QSO, bipolar beaming, 
errors in the continuum fits to the background galaxies' 
spectra, peculiar velocities, gradients in the matter den- 
sity, and so on) should also be symmetric about the QSO, 
at least in cosmic average, and so they will cancel out in 
a front-to-back comparison of the HI absorption. No sys- 
tematic that I can think of will have a shape that resembles 
the time delay surface in much detail. 

The closest candidate might be the gradual expansion 
of the universe, which causes the intergalactic material be- 
hind the QSO to be slightly denser than the material in 
front, changing the recombination time and producing a 
slight systematic gradient in the intergalactic HI density 
along the line of sight. This gradient could mimic a bright- 
ening QSO. The change in expansion scale factor as light 
traverses a region of proper length 40 Mpc at z = 3.0 is 
- 4% (fi M = 0.3, n A = 0.7, h = 0.65), and the cor- 
responding change in mean transmissivity of the IGM is 
A/ ~ 0.035 (e.g., McDonald et al. 2000). Af - 0.035 
is comparable to the smallest radiation-related change we 
might hope to measure (see below) , so this effect cannot be 

4 It may help to illustrate this point with a concrete example. 
Here is a crude recipe for reducing the noise in the control bins: fit 
a low-order polynomial to a plot of the mean transmissivity in each 
control bin versus the bin's distance to the QSO, then use the value 
of this function in each bin as the control transmissivity, rather than 
the measured transmissivity itself. Since we know a priori that the 
mean transmissivity in the control bins ought to be a smoothly and 
monotonically declining function of distance to the QSO, deviations 
from this behavior must be noise. They can be largely removed with 
the polynomial fit. Maximum-likelihood fitting of a surface to the 
unbinned data is a more sophisticated implementation of the same 
idea. 



ignored. One way to shrink it to insignificance is to divide 
each measured transmissivity by (f){z), the global mean 
transmissivity at its redshift. I will assume below that 
every transmissivity / has been rescaled in_a way (e.g., 
divided by (f)(z) and then multiplied by (/)(3) = 0.67) 
that makes this effect negligible. 

5. UNCERTAINTIES 

Suppose that we have divided the intergalactic volume 
surrounding the QSO into different spatial bins, each one 
illuminated by the light emitted by the QSO at a different 
period ti < t < ti + At in its past, and that we would 
like to use the mean Lyman-a transmissivity / = (e~ rLya ) 
in each bin to measure how the QSO's ionizing luminos- 
ity has evolved over time. How reliably can we do this? 
§ 3 discussed the relationship between / and the intensity 
of the radiation field. This section is concerned the un- 
certainty in our estimate of /. The interpretation of / is 
subject to further uncertainties, primarily related to our 
ignorance of the QSO's true proper distance z and to our 
assumption that its ionizing radiation is emitted isotrop- 
ically. The discussion of these will be deferred until §§8 
and 9. 

5.1. Arithmetic 

Consider a single spatial bin, for example the volume 
marked A in figure 3, that is being photoionized at the 
unknown rate T^^a- Let / Q b s be the mean Lyman-a trans- 
missivity we measure along the sight-line segments that 
pass through A, and let /true be the mean transmissivity 
that would be measured in an arbitrarily large intergalac- 
tic volume subjected to the same radiation field. / b s and 
/true will differ because (1) our spectra of the background 
sources are noisy, (2) we cannot measure the transmissiv- 
ity throughout the volume A but must instead rely on an 
uncertain guess from the few thinly scattered probes that 
background sources supply, and (3) A is not necessarily a 
fair sample of the universe and could have a mean trans- 
missivity that differs from /true f° r reasons unrelated to 
the intensity of ionizing radiation, e.g., if the random fluc- 
tuations from inflation gave it an unusually high or low 
density of matter. 

For simplicity I will obtain estimates of the variance 
from (1)— (3) by approximating the actual spatial bin A 
as the volume V that lies between two identically shaped 
parallel surfaces that are separated by distance l z (see fig- 
ure 4). Generalizing these results to the true geometry is 
conceptually trivial. In the simplified case, the variance 
due to (1) is simply 



ti 



fti 



(9) 



if we have spectra of identical signal-to-noise ratio S for 
each of the N background sources. 5 The variance due to 
(2) is related, in a way that is specified below, to the vari- 
ance of the mean transmissivity among randomly placed 
sightline segments of length l z , 



d 3 kP(k) 



sm(k z l z /2) 



k z l z /2 



(10) 



5 cS should be calculated for wavelength bins whose size corre- 
sponds to the depth l z of volume V, of course. 
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and similarly the variance due to (3) is related to the vari- 
ance of the mean transmissivity among randomly placed 
volumes of shape V, 



4 



-(h)' J 



d 3 rP(k)\W s (k x ,k v ) 



sin(M z /2) 



k z l z /2 



(11) 

Here -P(k) is the power-spectrum of transmissivity fluctu- 
ations and Ws(k x , k y ) is the Fourier transform of surface 
S. Equations 10 and 11 are both special cases of Parseval's 
relationship 

a2= (i)V rf3fcP ( k )l^( k )! 2 ( 12 ) 

between the power-spectrum and the variance of a random 
field that has been averaged over a volume whose shape 
has the Fourier transform IT(k). 



N background sources with 
positions specified by C(x,y) 




surf ac 



to earth 

Fig. 4. — Nomenclature for § 5. The intcrgalactic volume V is 
pierced by lincs-of-sight to many background sources. V is a general- 
ized right cylinder of height I z whose cross-section has the arbitrary 
shape S. Light from the background sources at the top of the figure 
moves parallel to the z axis, and normal to S, as it makes its way 
to towards earth, far below the bottom of the figure. The N back- 
ground sources have (x, y) positions that are randomly distributed 
within S. 

The first source of uncertainty, measurement errors, will 
simply add in quadrature to the others. To understand 
how the second and third sources contribute to the to- 
tal uncertainty, consider a given set C of random galaxy 
positions which can be represented as a sum of 

Dirac delta functions, C = (l/N) J2j ${ x ~ x j)^(v ~ Vj)- 6 
(The factor of 1 /N indicates our decision to average rather 
than sum the mean transmissivities from different galax- 
ies.) Equation 12 implies that the variance of / (i.e., of the 

6 Since galaxies are spatially extended, C should actually be the 
sum of functions with finite width, not the sum of delta functions. 
The delta-function approximation assumes that the bright part of 
each galaxy is small compared to the coherence length of intcrgalactic 
absorption. This appears to be the case at z ~ 3. 



transmissivity averaged over each of the N sightline seg- 
ments of length l z that pass through volume V) for this 
arrangement of background sources is 



a2 = w?l 



d 3 kP(k)\C k {k X ,ky)\' 



sm{k z l z /2) 



k z l z /2 



where 



|CW = ^E ex P[*Mr*-r m )] 



(13) 



(14) 



Lm 



is the powerspectrum of C, r rn = x m x + y m y is a vec- 
tor specifying the x,y position of the mth background 
source, and x and y are the usual unit vectors. Different 
arrangements of background sources will result in some- 
what different variances, but the average variance among 
all sets C of random galaxy positions is given by equa- 
tion 13 with |C fc (k)| 2 replaced by the expectation value 
(|C fc (k)| 2 ). The latter can be calculated by integrating 
over all random galaxy positions that lie behind the vol- 
ume V: 



<|C fe (k)| 2 } 



1 



E 



dxi dyi dx m dy m e lk (r; *" m) 



dxi dyi dx m dy m (15) 




The N(N — 1) terms in this sum with I ^ m are each 
equal to \Ws(k x , k y )\ 2 7 the powerspectrum of 5, while 
the N terms with I = m are each equal to unity. Sub- 
stitution into equation 13 shows that the expected vari- 
ance from the second and third sources of uncertainty 
is a linear combination of the variance in a volume of 
shape V and the variance along a line segment of length 
a 23 = ^W~ a v + W^f- T ms pleasantly simple re- 
sult makes intuitive sense. As N — > oo, sightlines to 
the background objects sample almost every part of the 
volume V. The mean transmissivity along the sightlines 
approaches the mean transmissivity within V arbitrarily 
closely and / bs becomes as reliable an estimator of / true 
as the volume-averaged transmissivity itself. In the op- 
posite limit, N = 1, the volume is pierced by a single 
skewer. In this case the footprint S of the volume (and 
therefore ay) becomes irrelevant; the mean transmissivity 
along the skewer is determined by physics on the other 
side of the universe, not by anything that influences S 
(e.g., our choice of telescope, instrument, pointing, etc.), 
and there is no reason that the variance of the transmis- 
sivity along this sightline segment should be any different 
than the variance of a similar segment that is randomly 
placed. 

Adding in the variance from measurement errors to the 
above expression for a\ 3 , we arrive at our final expression 
for the total variance of / : 



4 = 



N — 1 



1 



-oy 



Jtl 



(16) 



It remains to find numerical values for the constants of 
and ay. 

7 in the limit of weak angular clustering 
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High signal-to-noise QSO spectra can provide a robust 
empirical estimate of the variance of the mean transmis- 
sivity along a line segment of length l z . Figure 5 shows the 
value I find, as a function of l z , for the primary sample of 
7 QSOs at redshift z <~ 3 discussed in Adelberger et al. 
(2003). 




1 10 

I / h -1 comoving Mpc 



special cases of equation 11, 



ov = 



J dk R dk z k R P(k R , fc z ){ [ 



sm(k z l z /2) 



k z l z j2 



2Ji(fe fl i2)-|2 



k R R 



}} 



(19) 



or 



■ an 



dk R dk z k R P(k R ,k z ){ 
2R J 1 (k R R ) 



sm(k z l z /2)\- 



k z l z /2 
2R i J l {k R R^ 



k R {Rl-R}) 



] }, (20) 



will produce an estimate of cr v /af that is appropriate 
when V is a cylinder of radius R and height I (equation 19) 
or a cylindrical annulus of inner radius i£j, outer radius R , 
and height I (equation 20). Since the time-delay surface 
(§ 4) has rotational symmetry about the z axis, cylindrical 
bins are a natural choice for the analysis. 



Fig. 5. — The r.m.s. variation of the normalized transmissivity 
f/{f)(z) along skewers of different lengths. Points with error bars 
arc for the 7 primary QSOs from Adelberger et al. (2003); the shaded 
curve is the prediction of equations 17 and 10. The observed values 
at the largest I are unreliable due to uncertainties in the continuum 
fitting. 

ay is more difficult. It depends on the three-dimensional 
power-spectrum of the Lyman-a forest, which has not been 
measured due to a lack of close QSO pairs. A rough esti- 
mate of Oy can be obtained by scaling from of , a v = flaf , 
with P a constant that can be estimated by assuming a 
shape for the powerspectrum and numerically integrat- 
ing equations 10 and 11. According to McDonald (2003), 
the transmissivity powerspectra found in numerical sim- 
ulations the high-redshift intergalactic medium have the 
form 




1 10 
k z x h -1 comoving Mpc 



Fig. 6. — The one-dimensional power-spectrum of the Lyman-a 
forest. Points are from the observations of McDonald et al. (2001); 
the shaded curve is the prediction of equation 17. 



P(k R , k z ) oc [l + (3kl/k 2 \ P L (k)D(k R , k z ) (17) 

where k z and k R specify the wavenumber in polar coordi- 
nates, k 2 = k z + k R , Pl(/c) is the linear powerspectrum of 
cold dark matter, 



|Y k \"m / k \ a v ( k z \«- 
D{k R ,k.) -{-) -(— ^) 



(18) 



and k v (k) = (1 + k/k' v ) a -k v0 . McDonald (2003) sup- 
plies values for the numerical constants in these equa- 
tions appropriate to redshift z ~ 2. Adjusting his values 
slightly, to /? = 1.73, k ni = 5.2, a nl = 0.673, k p = 8.51, 
a p = 1.14, fc^o = 0.603, k' v = 0.085, a' v = 0.432, a v = 1.59, 
and adopting a T = 0.2 linear powerspectrum with a 
Pb{k) oc k 95 long-wavelength limit (Bardeen et al. 1986), 
I find a transmissivity powerspectrum that correctly pre- 
dicts both the observed dependence of of on l z at redshift 
z = 3 (figure 5) and the observed one-dimensional trans- 
missivity powerspectrum P 1D (k z ) = 2ir f dk R k R P(k R , k z ) 
at redshift z = 3 (figure 6). This model powerspectrum is 
presumably a reasonable approximation to the true pow- 
erspectrum. Integrating it numerically in the following 



A more sophisticated treatment would take into account 
the changes to the power-spectrum that will accompany 
changes in the ionizing radiation intensity. That will likely 
require numerical simulations and is beyond the scope 
of this paper. However, since changes in the radiation 
field alter the neutral fraction but not the temperature 
of intergalactic gas (see § 2), they will have a stronger 
effect on the amplitude of the powerspectrum than on 
its shape. Amplitude-independent results (e.g., our es- 
timate of & v /<jf) may not be disastrously affected. Fig- 
ure 7 presents some evidence in favor of this assertion. 
I converted published Lyman-a Voigt profile lists for the 
QSOs HS1946+7658 (Kirkman & Tytler 1997, z = 2.994), 
Q0636+680, Q0956+122, Q0302-003, and Q0014+813 (Hu 
et al. 1995, z = 3.180, 3.288, 3.294, 3.366) into a list of 
each QSO's Lyman-a optical depth as a function of red- 
shift, divided all the optical depths by a constant b to 
mimic a change in the ionizing radiation intensity, then 
calculated the dependence of oi on I for various values of 
b. The top panel of figure 7 shows the result. The am- 
plitude of (7i changes significantly with 6, but its shape, 
which is sensitive to the shape of the powerspectrum, does 
not. This can be seen more clearly in the bottom panel, 
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where the change in amplitude has been crudely removed 
by scaling each curve according to the relationship 

a t a [1 - /(6)] - 43 (21) 

where f(b) is the mean transmissivity after dividing the 
actual optical depths by b. Once this scaling is removed, 
the g i curves have nearly the same shape over the range 
1 £ I ^ lOh^ 1 comoving Mpc that is most important in 
our calculations. 
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Fig. 7. — Top panel: <Ji(l z ) (the r.m.s. transmissivity fluctua- 
tion along a line segment of length l z ) recalculated after dividing 
the observed Lyman-a optical depths in the QSOs HS1946+7658, 
Q0636+680, Q0956+122, Q0302-003, and Q0014+813 by the num- 
ber b written next to each curve. Bottom panel: same as top, except 
the dependence of the amplitude on the mean transmissivity f(b) 
has been crudely removed by dividing each curve by (1 — f(b))° 3 . 
The differences in shape revealed by this panel are small over the 
distance range 1 ^3 I lO/i^ 1 comoving Mpc that is most important 
to our analysis. 

5.2. Commentary 

It is worth considering the relative sizes of the terms in 
equation 16 before we move on. Suppose for concreteness 
that we are interested in measuring changes in the QSO's 
luminosity over lOMyr time-scales, so we have estimated 
the mean transmissivity in bins of depth l z ~ 1.5 proper 
Mpc, and suppose further that the radius R of our bins 
was chosen to fill the 40' field-of-view of a mosaicked CCD 
camera, R <~ 10 proper Mpc at z — 3 for = 0.3, 
Q\ = 0.7, h = 0.65. In comoving units the bin has R = 
2Qh~ 1 Mpc and l z = 3.9/i _1 Mpc, which implies ct; ~ 0.17 
(figure 5, for / = 0.67), crf/a v ~ 20, and consequently 
oy ~ 0.04. Inserting these numbers into equation 16 leads 
to the two primary conclusions of this section: 

(A) The spectra of the background objects do not need 
to be very good. The second and third terms in equa- 
tion 16 contribute equally to the total uncertainty when 
the signal-to-noise ratio in a spectral segment of length 



l z ~ 3.9/1" 1 comoving Mpc (~ 7A) is S = f/ai ~ 4. Ob- 
taining better spectra cannot change the total uncertainty 
by much. Random variations in the intergalactic matter 
density make the mean HI absorption /i along a single 
sightline segment a poor indicator of the ionizing back- 
ground J v and consequently there is no need to measure 
/i with exquisite precision. 

(B) The ultimate limit on our uncertainty is set by cry 
and it is a limit that we will reach very quickly. For the 
example considered here, the first and second terms in 
equation 16 have the same size for N = 20. Increasing 
N further can reduce the second and third terms arbi- 
trarily but not the total uncertainty. In practice one will 
want to obtain samples many times larger than N = 20 
to help address the possibility that the QSO's radiation 
is beamed. As a result oy will probably be the only sig- 
nificant contributor to <jf in realistic cases. The bottom 
panel of figure 2 shows that cry is large enough to prevent 
us from detecting all but the coarsest changes in a QSO's 
luminosity. We can do nothing about this. Only a small 
region of the universe is bathed in the light that the QSO 
emitted at one period in its history; the mean density in 
this region will stray from the global mean to the extent 
that inflation requires; the intensity of the QSO's ionizing 
radiation cannot be measured if it influences the region's 
mean HI absorption by less. 

Although the numbers quoted above depend on our ar- 
bitrary choice for the bin size, neither oi nor ay is a strong 
enough function of R or I for the qualitative conclusions to 
change significantly as the bin size varies across its useful 
range. The change of variance with mean transmissivity 
(equation 21) might seem more important. If the QSO 
were bright enough to drive the mean transmissivity in a 
bin to / ~ 0.9, for example, rather than the value / <~ 0.67 
assumed in the preceding three paragraphs, 07 would de- 
crease to 0.10, cry would decrease to 0.02, and one would 
find that spectra of slightly higher S:N were desirable. In 
practice, however, these high transmissivities are unlikely 
to be reached anywhere except very near the QSO, and 
here the decrease in 07 and ay is largely offset by the 
small bin sizes that are required at small radii. This can 
be seen in the worked example of § 7, below. 

One qualification should be added to my claim (A) that 
it is a waste of time to obtain high S:N spectra. That is 
necessarily true_only if one is committed to using the mean 
transmissivity / as a probe of the radiation density J v . At 
sufficiently high signal-to-noise ratios, however, other op- 
tions are available. Consider an arbitrary function g of the 
line-of-sight HI density that has variance due to random 
fluctuations when the brightness b = J v j 'Jb g of the ioniz- 
ing radiation field is kept constant. The random fluctua- 
tions in g will prevent us from detecting relative changes 
in b of order Aln6 ~ \dlnb/dlng\a g /g. If g is the mean 
transmissivity / on a 3.9/i _1 comoving Mpc line segment, 
which is the only possibility we have treated so far, then 
dlnb/dlng ~ 5 (near g — 0.67) and <J g jg ~ 0.26, so the 
minimum detectable fluctuation in b will have Aln6 ~ 1.3. 
At high enough S:N it would be possible abandon / and 
adopt (say) the total HI column-density M along the same 
line segment as the probe of b. In this case, with g = J\f, 
we have dlnb/dlng ~ —1 and u g jg ~ 6, where <J g /g is 
the value that obtains among the 5 QSOs with z ~ 3 
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discussed at the end of § 5.1, so the minimum detectable 
fluctuation in b has Aln6 ~ 6. TV is evidently far inferior 
to / as an estimator of b. A better choice, letting g equal 
n, the number of detected Lyman-a forest lines on the 
segment, was advocated by Bajtlik, Duncan, & Ostrikcr 
(1988). If the column density distribution has the form 
P(iV ffl ) oc TYjjf, then dlnb/dlng ~ (1 - a) -1 . Among 
the same 5 QSOs, <r n /n ~ 0.5 for line segments of length 
3.9/i _1 comoving Mpc, so the minimum detectable fluctu- 
ation in b has Alnfe ~ 1 for the observed slope a ~ 1.5. 
Surprisingly n is only marginally better _than / as an esti- 
mator of J v . The mean transmissivity / along a segment 
of a low S:N galaxy spectrum can provide almost as strong 
a constraint on the intensity of the QSO's radiation as a 
high S:N spectrum analyzed with the standard approach 
of Bajtlik, Duncan, & Ostrikcr (1988). It would be inter- 
esting to extend this analysis from cr; to the more relevant 
quantity cry. Since the ratio oil ay depends on the pow- 
erspectrum, and different estimators g will have different 
powerspectra, it is not necessarily true that the best esti- 
mator for o\ will be best for cry. This method would be 
made far more powerful if one could finding an estimator 
that significantly reduces cry. 

6. FEASIBILITY 

The previous section glibly claimed that observational 
uncertainties can be made smaller than cosmic variance. 
This section considers the claim in more detail. The op- 
timal depth for the spatial bins is ~ 7 A, corresponding 
roughly to 10 Myr time resolution (see § 8), and for the 
fields-of-view considered here the cosmic variance will con- 
sequently be cry ~ 0.04 (see §§ 5 and 7). This is the level to 
which we must reduce the observational uncertainty a no i sc . 

6.1. Signal 

Reducing the random errors to the desired level cr no i so ~ 
0.04 is not much of a challenge. Figure 8 shows spectra 
with <~ 10A resolution of three galaxies in the field SSA22a 
(Steidel et al. 2003) that were observed for ~ 29000s with 
the blue LRIS spectrograph (Steidel et al. 2004) on the 
Keck I telescope. The redshifts and AB-magnitudes of the 
galaxies (z ~ 3.1 and G ~ 24.5, respectively) are shown 
on the plot. The spectra are preliminary reductions that 
were selected more-or-less at random from the sample of 
Shapley et al. (2004, in preparation). These spectra have 
signal-to-noise ratios in the Lyman-a forest, for 7A bins 
(~ 1.5 proper Mpc), of around 5-6. Averaging together 
^10 spectra of similar quality would reduce the random 
errors to the desired level. Roughly forty spectra would 
be required if the exposure time were 2 hours instead of 8. 

6.2. Continuum fitting 

Systematic errors in the continuum fitting are a source 
of greater concern. We are interested not in a spectrum's 
flux itself but rather in its implied Lyman-a transmissivity, 
which is the ratio between the observed flux and the con- 
tinuum flux, i.e., the flux that would have been observed 
in the absence of Lyman-a absorption. The uncertainty 
in a bin's mean transmissivity therefore has an additional 
term that arises from errors in the estimated continuum. 

The size of these errors depends on the method that is 
used to estimate the continuum level. Traditional meth- 
ods are poorly suited to the present case; they exploit the 
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Fig. 8. — Examples of the Lyman-a forest in the spectra of galax- 
ies at redshift z ~ 3. The data were taken with LRIS-B (Steidel et 
al. 2004) and were generously provided by A. E. Shapley. Redshifts 
and AB magnitudes are indicated; the names ("Cll" etc.) refer to 
the SSA22a catalog of Steidel et al. (2003). Asterisks mark the loca- 
tions of the interstellar absorption lines discussed in figure 10, below. 
Dashed lines indicated the wavelengths of Lyman-a and Lyman-/? at 
the redshift of the galaxies' interstellar absorption lines. 

occasional presence of spectral regions with little Lyman- 
a absorption, and these regions are rare and difficult to 
recognize in noisy, low-resolution galaxy spectra. Several 
authors (e.g., Hui et al. 2001) have presented alterna- 
tives that are more useful to us. Particularly simple is the 
method of Croft (2004) , in which the continuum level c(A) 
is estimated by simply smoothing each objects spectrum 
and scaling appropriately: c(A) ~ /so(A)//(2:), where 
/50(A) is the object's Lyman-a forest spectrum smoothed 
by a Gaussian with a = 50A and f(z) is the published 
mean transmissivity at redshift z = A/AL yct — 1, which has 
been calculated by other authors from more-sophisticated 
continuum fits to the observed spectra of numerous bright 
QSOs. 

Experimentation on the 7 QSOs with z <~ 3 in the sam- 
ple of Adelberger et al. (2003) shows that transmissivi- 
ties estimated with Croft's (2004) approach and the tra- 
ditional approach are very similar. Excluding the parts of 
the Lyman-a forest that fall on the QSOs' Lyman-a and 
Lyman-/3 emission lines, the correlation coefficient between 
the transmissivities estimated with the two approaches is 
r ~ 0.96-0.98. This implies that the r.m.s. difference be- 
tween them, (1 - r 2 ) 1/2 CTi, is ~ 0.03-0.05 for 7A bins with 
gi ~ 0.17 (see Figure 5). Since galaxies' and QSOs' con- 
tinua between Lyman-a and Lyman-/3 are similarly fea- 
tureless (sec, e.g., figures 8, 9, and 10), errors in continuum 
fitting should introduce a similar uncertainty in galaxies' 
Lyman-a forest transmissivities. Even if these errors were 
correlated from one galaxy to the next, and did not tend to 
cancel in the averaged transmissivity in each spatial bin, 
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their size would be no larger than the cosmic variance. In 
practice continuum errors from the Croft (2004) approach 
should cancel significantly, however. They arise, to a large 
extent, because the Lyman-a forest is not completely uni- 
form on ~ 100A scales and any deviations from uniformity 
are incorrectly interpreted as features of the continuum. 
The resulting r.m.s. continuum error averaged over a spa- 
tial bin can be calculated with the approach of § 5: 



^cont 



JV-1 ,2 ,1 



N 



N 



1/2 



(22) 



where a' v and <r' ; are given by equations 11 and 10 with 
sin(fc z Z z /2)/(fc z Z z /2) replaced by exp[— fc z cr z /2], where o z 
is the comoving distance corresponding to the 50A smooth- 
ing length. Since er z 3> l z , cr'v wm be smaller than <jy, 
u 1 1 will be smaller than of, and the uncertainty from this 
source of continuum errors will never dominate the total 
uncertainty. This is true even in the low signal-to-noisc 
limit, since in this case the uncertainty in the smoothed 
continuum will be dwarfed by the uncertainty in a single 
~ 7A bin. 
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Fig. 10. — Mean absorption vs. wavelength for galaxies at rcdshift 
z ~ 3. The shaded curve shows the mean spectrum of 811 galaxies at 
z ~ 3; the data are taken from Shapley et al. (2003). Lyman-a forest 
fluctuations cancel out in this average, producing a smooth shelf 
with mean transmissivity / ~ 0.67 between Lyman-a and Lyman- 
(i. Absorption lines in this shelf are produced by material intrinsic 
to the galaxies. This interstellar absorption can be safely ignored 
when its rest- frame equivalent width is less than ~ 0.3A(see text). 
Stronger lines have been marked with vertical labels naming the ion 
responsible for the absorption. Data at these wavelengths should 
be excluded when calculating the mean transmissivity in different 
spatial bins. Note that the resolution of this spectrum (~ 12A) is 
somewhat lower than the 7A resolution advocated in the text. 
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Fig. 9. — Comparison of galaxy and QSO continua in the Lyman-a 
forest region. The top panel shows the gravitationally lensed galaxy 
MS1512-cB58 (Pettini et al. 2002); the bottom panel shows the 
QSO Q1623 KP78 (Adelberger et al. 2004, in preparation). Both 
spectra were taken with the same instrument, though the signal-to- 
noisc ratio for the bright QSO is significantly higher. Asterisks in 
the top panel mark the locations of the interstellar absorption lines 
discussed in figure 10, below. 



6.3. Interstellar absorption lines 

Although galaxies' continua appear to be mostly fea- 
tureless between Lyman-a and Lyman-/3, there are some 
important exceptions: at a handful of wavelengths the 
galaxies' interstellar absorption lines are too strong to be 
ignored. These wavelengths are apparent in figure 10, 



which shows average observed absorption in region be- 
tween Lyman-a and Lyman- (3 in a sample of 811 Lyman- 
break galaxies (Shapley et al. 2003). Since the r.m.s. vari- 
ation in Lyman-a forest transmissivity in 7A bins is 
o i ~ 0.17 (Figure 5), interstellar absorption lines will have 
a non-negligible effect on the estimated transmissivity if 
their rest-frame equivalent width exceeds 7Axcr;/(l + z) ~ 
0.3A. In the mean spectrum of Shapley et al. (2003) there 
are 7 such lines between Lyman-a and Lyman-/?. The ex- 
istence of these interstellar absorption lines will not have 
a disastrous effect on the analysis, since each background 
galaxy will lie at a slightly different redshift, but one might 
as well eliminate their effect completely by masking the 
relevant portions of the spectrum. This will reduce the 
effective number of background galaxies by too modest an 
amount to significantly alter the significance of any con- 
clusions. 

7. SYNTHESIS 

We can now estimate how easily we will be able to de- 
tect ionization gradients produced by changes in a QSO's 
luminosity. This section works through the details for a 
single case, an isotropically radiating QSO of (t = 0) mag- 
nitude mgi2 = 18 whose ionizing luminosity varied with 
time according to the curve L(t) shown in the top panel 
of figure 11. An Om = 0.3, f^A = 0.7, h = 0.65 cosmology 
with the uniform ionizing background of equation 2 will 
be assumed throughout this section. 

According to equation 6, the radiation intensity as a 
function of position in the observed frame (i.e., the radia- 
tion intensity that was present when the photons from the 
background sources passed through each point) is 

where ti(R, z) and r oq are given by equations 7 and 8. The 
left panel of figure 12 shows this function for r oq = 8.27 
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Fig. 11. — Top panel: the QSO's ionizing luminosity as a function 
of time. This curve, generated for the example of section 7, is not 
intended to be a realistic model of a QSO's output. It features were 
chosen to illustrate a number of points made in the text. Bottom 
panel: the recovered mean transmissivity in bins whose edges align 
with contours of the time-delay surface. Points mark the mean trans- 
missivity calculated from figure 12; error bars show the Icr range that 
would be observed if figure 12 included cosmic variance, the dom- 
inant noise source. The solid curve shows the mean transmissivity 
than would have been observed if the QSO's luminosity were con- 
stant, L(t) = L(0). The dotted curve shows the mean transmissivity 
in the primed "control" bins (figure 12); it differs from the solid line 
due to the increase in the QSO's luminosity from t = to t = —10, 
though this effect would presumably be removed by a maximum like- 
lihood fit of the data to a surface. Decreases of the luminosity to 
and subsequent increases back to L(0) can be detected with mod- 
erate significance if they happened within the interval — 50 t ^ 
Myr. Only extremely large increases in the QSO's luminosity can be 
detected at earlier times. 

proper Mpc, which is appropriate to the QSO described in 
the preceding paragraph. Since the ionization and recom- 
bination times (§ 2) are short compared to the time for 
significant changes in L(t), the neutral fraction will be in- 
versely proportional to J{R 1 z), and the mean correspond- 
ing transmissivity f(R, z) can be derived from the curve 
shown in the top panel of figure 2. The result is shown 
in the right panel of figure 12. The panel shows the mean 
transmissivity that would be observed if one averaged re- 
sults from many identical QSOs; the actual transmissivity 
surrounding a single QSO would have significant variations 
about this mean, due primarily to random fluctuations in 
the density of intergalactic matter. 

Changes in the QSO's ionizing luminosity could be de- 
tected in various ways, but for now I will assume that one 
is aiming to detect the changes by looking for differences in 
the mean transmissivity /bin among bins whose edges trace 
contours of the time delay surface ti(R, z). One set of such 
bins is shown in the left panel of figure 12. Also indicated 
are symmetrically distributed "control" bins. These bins 
have shapes identical to the others, but are located on the 
opposite side of the QSO, in a region that is illuminated 
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Fig. 12. — Left panel: variations in the intensity of ionizing 
radiation in the observed frame for a QSO at redshift 2 = 3 with 
"1912 = 18 that has the history of outbursts shown in figure 11. Also 
shown arc sample bins used in the analysis at larger lookback times. 
Right panel: variations in the expected mean transmissivity implied 
by the radiation intensity of the left panel. The mean transmissivity 
would hew to these values only if results from numerous identical 
QSOs were averaged. A more realistic single realization would show 
large fluctuations in the transmissivity due to random variations in 
the density of intergalactic matter. The labeled regions surrounding 
the QSO are sample bins for the analysis at smaller lookback times. 



by the QSO's radiation at short time delays > t > —10 
Myr. The mean transmissivity in these bins provides an 
indication of how the mean transmissivity for z > would 
vary with position if the QSO's luminosity were always 
equal to its observed (t = 0) value(§ 4). Note that the 
requirement t > — 10 Myr for the control bins limits the 
size of their partner bins at smaller radii. 

Our ability to detect changes in the QSO's luminosity 
will depend on the uncertainty in the binned transmissiv- 
ity /bin- For a reasonable number of background sources 
( few dozen) this uncertainty will be nearly equal to the 
cosmic variance cry (§ 5). As discussed in § 5, the size of cry 
can be roughly estimated by (1) approximating each bin 
as a cylinder of radius i?bin and depth ?bi n , (2) calculating 
r.m.s. transmissivity fluctuation along a skewer of length 
^bin by interpolating from figure 5 and scaling according to 
the local expected transmissivity f(R, z) with equation 21, 
and (3) multiplying by /3 1 / 2 , where f3(R hin ,l hin ) = er y /cr ; 2 
is a constant that can be determined for each bin by nu- 
merically integrating equations 10, 19, and 17. The bot- 
tom panel of figure 11 shows the expected mean transmis- 
sivity in each bin along with the estimated lcr uncertainty 
calculated in this way. 

The main conclusions from figure 11 are (1) that only 
large changes in the QSO's luminosity will leave an imprint 
on the IGM that is easily detectable with this approach, 
and (2) that the minimum size of a detectable luminosity 
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difference increases rapidly towards earlier times. These 
conclusions are easier to apprehend in figure 13, which 
shows the 10-Myr moving-average magnitude mgi2 a QSO 
would need to have had at various lookback times to for 
its sudden death (or revival) to leave a detectable fossil 
record in the spectra of background galaxies. Measuring 
changes in a QSO's luminosity will be difficult for t £ — 50 
Myr and nearly impossible for t £ — 100 Myr. At smaller 
time delays the prognosis is good. 



20 



does not depend sensitively on the details of this approxi- 
mation. 
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Fig. 13. — The lOMyr running-average luminosity required at 
various lookback times for the flux in the corresponding lOMyr bin to 
be No above the background level. Luminosity changes among even 
the brightest QSOs will be difficult to detect at moderate significance 
for t < - 100 Myr. 
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Fig. 14. — The mean transmissivity and its uncertainty in the 
bins shown in the right panel of figure 12. Bins with the lightest 
shading arc bathed in light emitted by the QSO at — 3.5Myr < t < 0; 
bins with medium shading by light emitted at —10 < t < — 3.5Myr; 
bins with darkest shading by —20 < t < — lOMyr. The height of 
the bins corresponds to the mean transmissivity ±1<t. Bin names 
in this figure correspond to those in figure 12. Bins with 1 < r < 
2Mpc arc unlabeled in figure 12 for clarity. At fixed time delay, / 
declines with radius due to the 1/r 2 decline in flux from the QSO. 
Comparing / radius by radius clearly reveals the decline in Lit) 
(figure 11) from -lOMyr < t < to -20 < t < -lOMyr. The 
increase from — 3.5Myr < t < to —10 < t < — 3.5Myr would not 
be easy to detect; a decrease during the same times would have been 
detected with high significance. 



As t — > 0, the shape of the time-delay contours becomes 
increasingly incompatible with the goal of having symmet- 
rical control bins. For —20 ^ t ^ the best approach may 
be to abandon the binning altogether and simply find the 
maximum likelihood fit of the data to an appropriate sur- 
face. To give some indication of the uncertainty in the ra- 
diative history that would result, however, I will continue 
a binned analysis with the bins shown in the right panel of 
figure 12. The outermost bins (B-BBB) enclose regions 
with time delays —20 < t < — 10 Myr, the middle bins 
(A- AAA) enclose — 10 < t < —3.5 Myr, and the innermost 
(AB' etc.) enclose —3.5 < t < Myr. A complication for 
these small time delays is that a region with fixed t has a 
wide range of distances r to the QSO. Variations in / due 
to changes in 1 jr 2 could obscure the variations of interest 
from changes in the QSO luminosity L(t). My approach in 
this simplified analysis is to divide each time-delay region 
into a number of bins with similar values of r (e.g., bins 
B, BB, and BBB for -20 < t < -10 Myr). Figure 14 
shows the mean transmissivity within each of these bins. 
The uncertainty in the bins was calculated by approximat- 
ing them as cylinders (B, AB' , AABB', AAABBB') or 
cylindrical annuli (A- AAA, B-BBB) with the same vol- 
ume and roughly the same shape; the expected variance 



With these bins, the decrease in L(t) from —10 < t < 
Myr to —20 < t < — 10 Myr is detected with some signifi- 
cance at each radius and with high significance when the 
results from different radii are combined. If the decrease 
had happened at an earlier time — 10 < t < — 3.5 Myr 
it would have been detected with higher significance still. 
Instead the luminosity increased from — 3.5Myr < t < 
to —10 < t ;$ — 3.5. As the plot shows, increases in the 
luminosity at small look-back times are much harder to 
detect than decreases, since the sensitivity of / to Jq falls 
as / — ► 1 (figure 2). Fortunately decreases must be far 
more likely than increases: QSOs with m 9 i 2 = 18 lie on 
the steep bright-end of the luminosity function. 

8. LIMITS TO THE TIME RESOLUTION 

In the previous sections the data were placed into spa- 
tial bins whose depth l z = 1.5 Mpc gave us sensitivity to 
luminosity fluctuations time-scales of 10 Myr or greater. 
Ideally one would be able to detect fluctuations on any 
time scale. Could we have achieved significantly better 
time resolution by placing the data in bins with smaller 
l z ? The answer is no; this section explains why. 

If cosmic variance were the only problem, the bin depth 
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l z could be made arbitrarily small, ay (equation 19) is 
almost independent of l z for l z <C R, the case of interest, 
and so a bin that is infinitcsimally thin will have nearly 
the same cosmic variance as the adopted bins with l z = 1.5 
proper Mpc. 8 

Unfortunately the ability to obtain reasonably accurate 
estimates of the mean transmissivity in infinitcsimally thin 
bins is not the same as the ability to measure changes in 
the QSO's luminosity that happened on arbitrarily short 
time scales. Our estimate of a gas element's longitudi- 
nal separation z from the QSO will be inaccurate for two 
reasons: (1) we will not know the precise redshift of the 
QSO, and (2) our z positions are derived from redshifts 
and will be distorted from their true values by peculiar 
velocities. As a result the time delay to the element is 
uncertain. When the time delays to different elements are 
uncertain, we cannot combine elements with exactly the 
same delays into one bin; the various elements that make 
up a single bin will inevitably have a range of time de- 
lays. The minimum range of time delays in a bin is what 
limits our time resolution. The remainder of the section 
considers this limit in more detail. 

8.1. Uncertainty in the QSO redshift 

If the QSO's redshift is measured from the CIV emission 
line, the uncertainty in its systemic recession velocity will 
be u v ~ 510km s -1 (Richards et al. 2002), which corre- 
sponds to a positional uncertainty of a z ~ 1.75/ig 5 proper 
Mpc at z = 3 for Qm = 0.3, Qa = 0.7. The uncertainty 
can be reduced to 0.9/igg 1 proper Mpc (i.e., 270kms _1 ) if 
Mgll is used instead (Richards et al. 2002), and to O.S/ig^ 1 
proper Mpc (i.e., SOkms" 1 ) if [OIII] is used (Vrtilek & 
Carleton 1985). Radio observations of molecular emission 
lines could presumably reduce o~ z even further, but this is 
unlikely to benefit us much. The time resolution achiev- 
able with a z — 0.3 proper Mpc is t <~ 2a z /c ~ 2 Myr, and 
other effects prevent us from obtaining a resolution even 
this coarse. 

8.2. Thermal motions 

A firm lower limit to the time resolution is set by 
the thermal motions of the 20000 K intergalactic gas. 
The intergalactic hydrogen at a particular true z posi- 
tion will have an rms range of apparent z position of 
g z = {kT/m H ) 1/2 /H ~ 0.04 proper Mpc for il M = 0.3, 
Oa = 0.7, h = 0.65. We will not be able to measure 
changes in the QSO's luminosity that happen much more 
rapidly than the corresponding time scale t ~ 2az/c <~ 0.3 
Myr. This is unlikely to be the limiting factor in the anal- 
ysis. 

8.3. Streaming towards the QSO 

The effect of larger-scale peculiar velocities is more se- 
vere. First there is the average streaming motion towards 
the QSO, which can be crudely estimated as follows. If the 
scale dependence of QSOs' bias b is weak, then the mean 

8 The reason is that a cylinder with l z <C R in real space will 
have l z 2> R in Fourier space, and since the power P(k) is concen- 
trated near k ~ the variance is dominated by contributions from 
wavenumbers k near the origin. There is little power at the distant 
ends of the cylinder with large fc z , and altering the limits of the k z 
integral in equation 19 (i.e., changing the thickness l z of the cylinder 
in real space) does not affect the variance by much. 



matter overdensity at a distance r from a QSO will be 
roughly 8(r) = S,Q{r)/b, where £q is the correlation func- 
tion of QSOs. According to Croom et al. (2002), a corre- 
lation function of the form £g(r) = (r/r )~ 7 , r ~ 8.4/i _1 
comoving Mpc, 7 ~ 1.56 is appropriate for the brightest 
QSOs at any redshift. The variance of QSO number den- 
sity in cells of radius r ce u = 8/i _1 comoving Mpc is there- 
fore 0% = 72(r /r ce ii)'7[(3 - 7) (4 - 7)(6 - 7)2 7 ] * 1.3, 
(Peebles 1980 eq. 59.3), which implies a QSO bias at 
z = 3 of b ~ 4.5 if CIm = 0.3, Oa = 0.7, and the rms 
linear matter-density fluctuation in r = 8/i _1 Mpc spheres 
at z — is cr 8 = 0.9. Integrating over the correlation 
function shows that the mean matter overdensity within a 
comoving radius r c , S(r c ), is roughly 



6(r) 



3r7 7 



(3 - 7)6r ' 



(24) 



Now the proper radius r p of a spherically symmetric region 
with mean interior overdensity 8 evolves according to 



H 



dS 



3(1 + 5) dt 



(25) 



where H{t) is the Hubble parameter. This follows from 
the fact that concentric shells of matter do not cross un- 
til just before final collapse. Adopting the spherical Zel- 
dovich approximation for simplicity, the linear overdensity 
density 8 l will be related to the true overdensity through 
1 + 8 ~ (1 — <5l/3)~ 3 , which reduces equation 25 to 



H 



l + f-f(l + 8(r p )) 



1/3 



and shows that 



Az = f 



(R 2 + Z 2 ) 1 / 2 



1 - (1 + S(r p )) 



1/3' 



(26) 



(27) 



is the rough proper distance between the a volume ele- 
ment's true position (R, z) and the position we that we 
erroneously infer from assuming that it is at rest with re- 
spect to the Hubble flow. Here / = d\nD/d\na ~ Cl^(z) 
where D is the linear-growth factor and a is the scale fac- 
tor of the universe. Figure 15 shows 8 and Az as a function 
of distance. 

By itself the net streaming towards the QSO is not a 
major problem, at least for r ^ 5/i _1 comoving Mpc. Its 
primary effect is to make absorbing gas appear to be closer 
to the QSO than it actually is. Although this produces a 
slight systematic error in the lookback time assigned to 
each volume element, the error can be corrected to a large 
degree with simple formulae such as equation 27, and in 
any case slight inaccuracies in the times assigned to the x 
axis of figure 11 would not diminish its scientific value by 
much. 

8.4. Random peculiar velocities 

More troubling are random deviations around the net 
streaming motion. As a result of them, the gas that is 
illuminated by the QSO's luminosity at time t will lie on 
a complicated surface that wanders randomly around the 
parabolic time-delay surface shown in figure 3. Since par- 
ticles maintain their linear velocities long after the density 
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where 9 is the angle between r and the z axis, 
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This is a special case of a result derived by Gorski (1988). 
If the wavenumbers in the integrals are comoving, as is the 
convention, the rms error in comoving z separation due 
to peculiar velocities is <r z (r) = a Vz (r)/H. Inserting the 
T = 0.2 linear powerspectrum of Bardcen et al. (1986) into 
these equations, normalizing to er 8 = 0.9 at redshift z = 
(i.e., to erg — 0.29 at z = 3, appropriate for Qm = 0.3, 
Oa = 0.7), and integrating numerically, I find the values 
of a z shown in figure 16. 



Fig. 15. — The mean interior overdensity (solid line) and ap- 
proximate infall velocity (dashed line) as a function of distance 
from a bright QSO at z = 3. The infall velocity is expressed as 
Ar/r = (r ap parcnt - rtrue)/Vtrue, the relative shift in an object's 
position that the infall velocity would produce if the velocity were 
directed exactly towards the observer. The infall velocity is esti- 
mated with a simple version of the Zcldovich approximation that 
becomes inaccurate for 8 2> 1; sec text. 



field itself has left the linear regime, and since most of in- 
tergalactic space should be occupied by low density (i.e., 
uncollapsed) gas that is not far from the linear regime, lin- 
ear perturbation theory should provide a rough estimate 
of the typical size of these excursions. Let a Vz (r) be the 
rms difference in the z comoving peculiar velocities of two 
points separated by the vector r. In the linear regime the 
Fourier transform v(k) of the comoving peculiar- velocity 
field v(r) is related to the Fourier transform <5(k) of the co- 
moving density field <5(r) through v(k) = —iH fS(\s.)]s./k 2 , 
(e.g., Peebles 1980 equation 27.22) where H is the Hubble 
parameter and / ~ iVj^(z). Convolving v z (r) by the sum 
of a positive delta function at r'/2 and a negative delta 
function at — r'/2 produces a new random field whose value 
at each point is equal to the ^-velocity difference between 
the points r'/2 and — r'/2. The variance of this field is 
equal to a 2 z (r) , which can therefore be written, according 
to equation 12, as 

< (r) = ^ 3 4 ft(k)sin2 (^)- (28) 

Expressing a 2 ^ (r) in terms of one-dimensional integrals by 
converting to spherical coordinates and integrating over 
solid angle, one finds 

<(r) = 2[V>||(0) -^(r)cos 2 (0) - ^(r) sin 2 (0)], (29) 
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FlG. 16. — The peculiar velocity field. Peculiar velocities add 
random offsets to the apparent z positions of intcrgalactic gas and 
systematically shift them towards the QSO. Black points mark the 
true spatial positions of selected intcrgalactic volume elements; error 
bars show the expected lcr range of each volume element's apparent 
position. The mean offsets were calculated with equation 27 and the 
sizes of the ranges with equation 29. Contours in the background 
show the time delay to each position, starting with 1, 3, 10, 20 Myr 
and increasing in steps of 20 Myr thereafter. Because of peculiar 
velocities, gas with a range of time delays is mixed together in any 
spatial bin that is defined in the observed (redshift-space) frame. 
This limits the time resolution of the method. For time delays greater 
than ~ 20 Myr, the time resolution is limited to £ 10 Myr. At very 
small time delays the resolution is far better. Peculiar velocities are 
unlikely to make material with a delay of ~ 1 Myr appear to have a 
delay of 3 Myr, for example, so we should be able to use the method 
to detect QSO lifetimes of order 1 Myr even though we will not be 
able to distinguish a life time of 30 Myr from one of 31 Myr. 



8.5. Upshot 
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Random peculiar velocities are likely to be the dominant 
source of uncertainty in a volume element's distance to the 
QSO. The uncertainty in the element's proper z position 
is a z ~ 1 proper Mpc, which corresponds to a time-scale of 
t ~ 2(T z /c ~ 7 Myr. It might be possible, with very high 
signal-to-noise spectra for a very large number of back- 
ground sources, to trace and correct these distortions to 
the time delay surface. A number of interesting applica- 
tions would then be possible. With current technology, 
however, we will have accept that the mean transmissivity 
in any bin we devise will be sensitive to the QSO's lumi- 
nosity over a ~ 10 Myr range of lookback times. Time 
resolution significantly better than this does not appear 
to be achievable. 9 

9. SUMMARY AND DISCUSSION 

This paper showed that changes over time in the lumi- 
nosity of a QSO at redshift z ~ 3 will produce ionization 
gradients in the IGM and alter the Lyman-a forest ab- 
sorption spectra of background galaxies in an observable 
way. Because the density of detectable galaxies (1Z ^ 25) 
at z <~ 3 is high, <~ larcmin -2 , their absorption spectra 
can provide a detailed view of the ionization gradients. If 
an isotropically radiating QSO has an AB magnitude at 
912A of m 9 i2 = 18, significant decreases in its luminosity 
at larger look-back times will be detectable if they happen 
1 t £ 50 Myr before the time of observation. The time 
limits expand for brighter QSOs and shrink for fainter. In- 
creases in the QSO's luminosity over this time period will 
be harder to detect than decreases, but since mgi2 = 18 
corresponds to the steep bright end of the QSO luminos- 
ity distribution, they must be much rarer. § 7 sketches 
out the method and presents the uncertainties for a simu- 
lated QSO with known radiative history L(t); the section 
is aimed at those who want more detail but are reluctant 
to read the entire paper. 

The method gives us sensitivity to changes in a QSO's 
luminosity over a useful range of times. Statistical ar- 
guments mentioned in the introduction show that QSOs 
must change their luminosities significantly on time-scales 
t £ 100 Myr. If these changes happen on time-scales 
t < 0.1-1 Myr, a handful of QSOs in large (SDSS-sized) 
samples will show major brightness changes from one 
decade or century to the next (Martini & Schneider 2003). 
The method I have described cannot detect luminosity 
changes that happen on time-scales so short (§ 8), but is 
sensitive changes throughout the rest of the allowed range 
(1 £ t ;$ 100 Myr). Taken together, the two methods will 
be able to pin down the typical QSO lifetime in a robust 
and direct way. 

A number of other authors (e.g., Crotts 1989; Dobrzycki 
& Bechtold 1991; Moller & Kjaergaard 1992; Fernandez- 
Soto et al. 1995; Liske & Williger 2001; Jakobsen et al. 
2003; Schirber, Miralda-Escude, & McDonald 2004; Croft 
2004) have attempted to measure QSO lifetimes with a 
similar approach. Their results were ambiguous. The rea- 
son is that they used the absorption lines in a single back- 
ground QSO (or, in some cases, a handful) to search for 
ionization gradients around the foreground QSO. When 

9 Exccpt at very early times, > t > —20 Myr, when the posi- 
tional uncertainties become increasingly aligned with the time-delay 
contours; see figure 16. 



the QSO pair had a small projected separation, the analy- 
ses were confused by the high densities and large peculiar 
velocities near the foreground QSO (see, e.g., figure 15); 
when the projected separation was large the foreground 
QSO's weak effect on the IGM could not stand out above 
the cosmic variance. This paper's method sidesteps these 
difficulties by using numerous faint galaxies rather than a 
small number of bright QSOs as the background sources. 
As shown in § 5, the signal-to-noise ratio of the background 
objects' spectra does not affect the final result by much, 
but the number of background sources does. Choosing 
galaxies as the background sources is therefore the sensible 
approach. With numerous background sources, the weak 
effect of a QSO on distant intergalactic matter can be de- 
tected with reasonable significance, the complicated region 
closest to the QSO can be ignored altogether, and the pe- 
culiar velocity and density gradients at slightly larger dis- 
tances can be compensated by comparing to the amount 
of Lyman-a absorption in the large "control region" that 
is illuminated by light emitted by the QSO at t ~ (§4). 
The approach I have described should therefore be a sig- 
nificant improvement over previous work. 

Observers who would like to apply this approach in prac- 
tice should be aware that the optimal proper size of the 
observed region can be very different from the naive guess 
c^max with i max the maximum time delay of interest. As 
imax — ► intergalactic absorption at radii R 3> ci max be- 
comes increasingly important for the analysis, as figures 3 
and 14 show. Two results derived in § 5 are also relevant. 
They are discussed more fully in § 5.2, but may be sum- 
marized as follows: (A) Galaxy spectra with S:N ~ 4 per 
7A bin will be sufficient for measuring the QSO's radia- 
tive history with the stated precision. Obtaining better 
spectra will not improve the result by much. (B) Cosmic 
variance places a fundamental limit on the accuracy of the 
method. Only a small part of the universe is bathed in 
the light that the QSO emitted at a particular moment 
in its history, and the HI content of this region will stray 
from the global average for reasons that have nothing to 
do with the luminosity of the QSO. The only detectable 
luminosity changes are those large enough to alter the in- 
tergalactic HI content by more than its intrinsic random 
fluctuations cry. A corollary is that the number of back- 
ground sources docs not need to be incredibly high; one 
needs only enough to measure the mean transmissivity of 
a region with a precision similar to cry, and this is eas- 
ily achieved with spectra of a few dozen galaxies (but see 
below) . 

I have neglected an important complication in the dis- 
cussion so far. It is possible that the QSO's ionizing ra- 
diation will not emerge isotropically but will instead be 
focused into a bipolar beam with opening angle a ~ 90° 
(Barthel 1989). If the beam were pointed towards earth 
there would be almost no effect on the analysis (see fig- 
ure 12), but in the typical case the intergalactic volume 
that is affected by the QSO's radiation will be only ~ 30% 
as large as I have assumed. This will increase the un- 
certainty in the results. Shrinking the radius of one of 
our idealized cylindrical bins (§§ 5 and 7) until it contains 
only 30% of its previous volume will increase its uncer- 
tainty due to cosmic variance, ay, by a factor of ~ 1.5. If 
cosmic variance is still the dominant source of noise in the 
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smaller bins, the error bars in the / vs. t curve will increase 
by a similar factor. Fortunately this is not enough to pre- 
vent us from detecting changes in the QSO's luminosity for 
— 50Myr t (figure 11). Moreover, even with the enlarged 
error bars it should be easy to distinguish intergalactic 
volumes that are illuminated by the QSO's ionizing radi- 
ation from neighboring regions with transmissivities close 
to the global mean of / ~ 0.67, at least at small radii 
where the QSO's radiation is most intense (see, e.g., fig- 
ure 14, which shows that the QSO's influence at small radii 
should be detected with high significance) . This statement 
is independent of the QSO's radiative history as long as it 
has been shining for more than ~ 10 Myr, since so large 
a volume is (potentially) illuminated by light emitted at 
— lOMyr < t < 0. If one were still worried about the possi- 
bility of beaming, QSOs could be chosen for study only if 
they have characteristics that suggest their beam is likely 
to be pointed towards the earth (e.g., BAL QSOs or radio- 
loud QSOs). I suspect that it will be more interesting to 
observe many types of QSOs and use this approach as a 
direct test of unified models of AGN. 

In any case, the possibility of beaming makes it clear 
that the ideal number of background sources is in fact 
several times larger than the arguments of point (B) sug- 
gest (see above). We would like to be cosmic- variance lim- 
ited in even the possibly small fraction of the field that is 
struck by the QSO's ionizing radiation. Several hundred 
background sources would be ideal. Fortunately multi- 
object spectrographs with the required large field-of-view 
(e.g., IMACS; Dressier et al. 2003) can obtain this many 
spectra with a small number of slitmasks. Achieving the 
necessary signal-to-noise ratio for the background galax- 
ies at z 3 is less daunting than it sounds, since the large 
field-of-view allows one to pick sources from the bright end 
of the luminosity distribution. 

I should mention in closing that the results of this pa- 
per could be extended or improved in a number of ways. 
The expected cosmic variance cry played a large role in 
the analysis, yet was estimated from an imperfect model 
of the three-dimensional transmissivity power-spectrum. 
Better numbers should be derived from numerical simu- 
lations. It would be interesting to know if measuring the 
radiation intensity with something other than the mean 
transmissivity / would let us achieve tighter constraints 
on the QSO's radiative history before we were limited by 
cosmic variance. I assumed that the QSO would have red- 
shift z = 3.0, but in fact observers could choose to observe 
QSOs at any redshift 2 & z ;$ 4 where the Lyman-a for- 
est is visible from the ground and reasonable numbers of 
background galaxies can be identified. Finding the QSO 
redshift that minimizes the required observing time would 
allow one to study a larger number of targets. This work 
is admittedly unfinished. I hope to have shown that it is 
worth pursuing further. 



George Becker and Luis Ho were a useful sources of 
information about AGN and the uncertainty in different 
QSO-redshift estimators. Esther Hu helped me find elec- 
tronic version of her Voigt-profile line lists. Chuck Steidel 
fielded my random questions like a pro. Alice Shapley's 
email responses to my questions were always prompt and 
enlightening, and in addition she graciously allowed me to 



show the data in Figure 8. Rob Simcoe had helpful advice 
on several subjects. Paul Martini offered useful comments 
on an earlier draft. Andrew Carnegie's generosity put the 
food on my table. It is a pleasure to acknowledge the hos- 
pitality of Las Campanas Observatory during a long stay 
in Dec 2003 when this paper was begun. 
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